skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Thatcher, D"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Greenland stadials and interstadials (GS/GI) were millennial climate oscillations during the last glacial period that were originally identified in Greenland ice cores but that have been correlated with environmental change around much of the globe, including in monsoon regimes, with enhanced monsoon rainfall coincident with North Atlantic warming. Hydroclimate variability associated with GS/GI have been investigated in detail using terrestrial (primarily oxygen isotopes in stalagmites) and marine records, particularly for the Southeast Asian monsoon. However, a considerably smaller number of terrestrial records preserve these events in the Indian summer monsoon (ISM), the primary water source for ~2 billion people across South Asia. Here we present the first glacial-age speleothem stable isotope time series from Nepal, located in the ISM regime. UK-1 is a 187 mm tall aragonite stalagmite from the Pokhara Valley of central Nepal, ~150 km west of Kathmandu. The chronology of UK-1, which was established by 8 U/Th dates, all of which fall in stratigraphic order (within the errors), reveals continuous growth from 34,350-31,500 yr BP (Marine Isotope Stage 3); age uncertainties average ±84 yr. Stable isotope samples were measured every 1 mm, corresponding to a temporal resolution of 18 yr. Oxygen isotope ratios range from -5.6‰ to -7.6‰, and share the same timing and structure as Greenland (inter)stadials GS/GI 6 and 5.2 in the NGRIP record. We interpret this as reflecting an amount effect response to a strengthened ISM driven by more (less) poleward migration of the intertropical convergence zone during periods of northern hemisphere warming (cooling). This clear millennial signal in UK-1 is a somewhat unexpected result given that amount effects in oxygen isotopes in precipitation are weak (R^2=0.1) in this area today. UK-1 carbon isotope ratios range from -3‰ to -6‰ (excluding a small number of negative spikes) and exhibit variability coarsely similar to the NGRIP record, with lower (higher) values generally corresponding to GI (GS), possibly due to prior calcite precipitation in voids above the cave concomitant with changes in precipitation. Some periods of antiphasing between carbon and oxygen are also apparent and may reflect flushing of soil carbon dioxide during particularly wet phases. 
    more » « less
    Free, publicly-accessible full text available April 1, 2026
  2. Free, publicly-accessible full text available December 1, 2025
  3. Over the late Holocene, a variety of hydroclimate-sensitive proxies have identified substantial, multidecadal changes in Indian summer monsoon (ISM) precipitation, the most prominent of which is the “4.2 ka event”. This interval, dated to ~4.2-3.9 ka, is associated with severe droughts across South Asia that are linked to societal change. Given the absence of the 4.2 ka event in polar records, the 4.2 ka event is generally associated with low latitude forcings, but no clear consensus on its origins has been reached. We investigated the ISM response to the 4.2 ka event through analysis of aragonite stalagmites from Siddha cave, formed in the lower Paleozoic Dhading dolomite in the Pokhara Valley of central Nepal (28.0˚N, 84.1˚E; ~850 m.a.s.l.). The climate of this region is dominated by small monthly variations in air temperature (21±5˚C) but strong precipitation seasonality associated with the ISM: ~80% of the annual 3900 mm of rainfall occurs between June and September. High uranium and low detrital thorium abundances in these stalagmites yield precise U/Th ages that all fall within stratigraphic order. These dates reveal continuous growth from 4.30-2.26 ka, interrupted only by a hiatus from 3.27-3.10 ka. Overlap with coeval aragonite stalagmites reveals generally consistent trends in carbon and oxygen isotope ratios, suggesting that these stalagmites reflect environmental variability and not secondary (e.g., kinetic) effects. Many stalagmite-based paleomonsoon reconstructions rely on oxygen isotope ratios, which track amount effects in regional rainfall. However, our on-going rainwater collection and analysis program, as well as a previous study conducted in Kathmandu, 120 km the east of Siddha cave, reveals that amount effects in precipitation are weak in this region, particularly during the monsoon season, and thus we rely instead on carbon isotope ratios, which have been demonstrated to track site-specific effective precipitation. Siddha cave stalagmite carbon isotopes, in contrast to other South Asian proxy records, indicate that ISM rainfall increased at Siddha cave from 4.13-3.91 ka. As a further test of this result, we analyzed uranium abundances in the section spanning 4.3-3.4 ka. Uranium serves as an indicator of prior aragonite precipitation and thus of hydroclimate, and like carbon isotopes, suggests increased ISM rainfall coincident with the 4.2 ka event. This precipitation anomaly is nearly identical in timing and structure but anti-phased with stalagmites from Mawmluh cave, northeastern India. We investigated the climatic origins of this precipitation dipole using observational data from the Global Precipitation Climatology Centre (GPCC) and Hadley Center Sea Ice and Sea Surface Temperature (HadISST) products. Preliminary spatial composites suggest that large precipitation differences between Mawmluh and Siddha caves are associated with SST anomalies in the equatorial Pacific. Additionally, superposed Epoch Analysis shows relatively rapid eastern Indian Ocean cooling during the summer monsoon season coeval with large precipitation differences between these sites. Our findings lend support to a tropical Indo-Pacific origin of the 4.2 ka event. 
    more » « less
  4. Parts of the northwest Atlantic Ocean, including the Gulf of Maine, along the northeastern coast of the United States, are warming at a rate as much as four times faster than the global ocean, according to instrumental and satellite records. In addition to the longer-term temperature increases, the frequency and severity of marine heat waves have been increasing. Combined, the ecological impacts are numerous and concerning, yet our understanding of past climate in this region is incomplete due to sparse and short-duration instrumental records. Here we present annually resolved oxygen isotope records from the marine bivalve, Arctica islandica, from five locations ranging from Jonesport and Seguin Island in the Gulf of Maine to Long Island, New York, Cape May, New Jersey, and Ocean City, Maryland in the Mid-Atlantic Bight, a span of over 870 km along the Atlantic coast. Several of the isotope records span the last 100 years or more and all records show coherent, substantial warming since at least 1980 CE. The level of warming indicated in the shell oxygen isotopes is comparable to the 0.5 °C per decade (1980-2020 CE) warming also shown in the instrumental record of sea surface temperature from Boothbay Harbor along the central coast in the Gulf of Maine. These five spatially distant isotope records span different oceanographic conditions and dynamics, including water mass sources, yet they all indicate a substantial warming in recent decades, likely related to increased anthropogenic warming. Beyond reconstructing seawater temperature prior to instrumental records, a major goal of this work is to disentangle the global warming signal from these records to better understand the underlying ocean dynamics also influencing these records. 
    more » « less
  5. Coastal systems can exhibit large variability in pH compared to open marine conditions, thus the impacts of ocean acidification (OA) on their resident calcifying organisms are potentially magnified. Further, our understanding of the natural baseline and variability of pH is spatially and temporally limited in coastal settings. In the few coastal locations that have been monitoring seawater pH, records are generally limited to <10 years and are thus unable to provide the full range of centennial to decadal natural variability. This is the case for the Gulf of Maine (northwestern Atlantic), a highly productive region of strategic importance to U.S. fisheries, that is facing multiple environmental stressors including rapid warming and threats from OA. Paleoceanographic proxy records are therefore much needed in this region to reconstruct past pH conditions beyond instrumental records. A clear candidate for this is the boron isotope (d11B) pH proxy provided the d11B sensitivity to pH in long-lived shallow water marine carbonates can be established. To this end, we grew juvenile and adult Arctica islandica (ocean quahog) in flowing seawater tanks for 20.5 weeks in controlled pH (7.4, 7.6, 7.8 or 8.0 (ambient) ± 0.02) and temperature (6, 9 or 12 ± 0.56 °C) conditions at Bowdoin College’s Schiller Coastal Studies Center, Harpswell, Maine (USA). The clams were stained twice with calcein and supplemented with food (Shellfish Diet) throughout the experiment to ensure suitable growth. New shell growth (average 67% increase in maximum shell height and 522% increase in buoyant weight across all treatments), constrained by calcein markings, were sampled for boron isotope analysis (d11B) to determine if shell d11B varied as a function of pH similar to many other calcifying organisms. The results of the culture experiment will yield whether or not Arctica islandica preserves seawater pH information in their shells. If so, the transfer function relating shell d11B to pH will be used to hindcast pH in the central coastal region of the Gulf of Maine during recent centuries. Alternatively, if the shell d11B signal is independent of ambient seawater pH, this may reveal the capacity of Arctica 
    more » « less
  6. Coastal systems can exhibit large variability in pH compared to open marine conditions, thus the impacts of ocean acidification (OA) on their resident calcifying organisms are potentially magnified. Further, our understanding of the natural baseline and variability of pH is spatially and temporally limited in coastal settings. In the few coastal locations that have been monitoring seawater pH, records are generally limited to <10 years and are thus unable to provide the full range of centennial to decadal natural variability. This is the case for the Gulf of Maine (northwestern Atlantic), a highly productive region of strategic importance to U.S. fisheries, that is facing multiple environmental stressors including rapid warming and threats from OA. Paleoceanographic proxy records are therefore much needed in this region to reconstruct past pH conditions beyond instrumental records. A clear candidate for this is the boron isotope (d11B) pH proxy provided the d11B sensitivity to pH in long-lived shallow water marine carbonates can be established. To this end, we grew juvenile and adult Arctica islandica (ocean quahog) in flowing seawater tanks for 20.5 weeks in controlled pH (7.4, 7.6, 7.8 or 8.0 (ambient) ± 0.02) and temperature (6, 9 or 12 ± 0.56 °C) conditions at Bowdoin College’s Schiller Coastal Studies Center, Harpswell, Maine (USA). The clams were stained twice with calcein and supplemented with food (Shellfish Diet) throughout the experiment to ensure suitable growth. New shell growth (average 67% increase in maximum shell height and 522% increase in buoyant weight across all treatments), constrained by calcein markings, were sampled for boron isotope analysis (d11B) to determine if shell d11B varied as a function of pH similar to many other calcifying organisms. The results of the culture experiment will yield whether or not Arctica islandica preserves seawater pH information in their shells. If so, the transfer function relating shell d11B to pH will be used to hindcast pH in the central coastal region of the Gulf of Maine during recent centuries. Alternatively, if the shell d11B signal is independent of ambient seawater pH, this may reveal the capacity of Arctica islandica to regulate internal calcifying fluid chemistry and their resilience to OA. 
    more » « less
  7. Arctica islandica (ocean quahog), a commercially-important, long-lived bivalve species, is abundant on much of the northeastern United States continental shelf. Several recent studies have noted increases in growth rates of these clams over the last 200 years at some locations in the southern Mid-Atlantic Bight region whereas growth rates at sites farther north have remained constant through time. It has been suggested that these changes in growth rate are related to warming in the more southerly sites. However, a direct comparison between site-specific bottom-water temperatures and A. islandica growth rates has not been done. We present oxygen isotope data measured in Arctica islandica shells, a proxy for seawater temperature, paired with simulated temperature from high-resolution ocean model output to investigate the relationship between A. islandica shell growth rate and bottom water temperatures throughout the northeastern United States continental shelf. The relationship between oxygen isotopes and growth rate in A. islandica is assessed at several locations, including the continental shelf offshore New Jersey and Long Island, and the Georges Bank region. Bottom water temperature trends at these locations are further assessed using the VIKING20X ocean model, which uses JRA55-do (55-year Japanese Atmospheric Reanalysis for driving ocean-sea-ice models) atmospheric forcing from 1958 to present and nests a 1/20° Atlantic Ocean in a 1 ⁄ 4° global domain. The results of this work have implications for the ocean quahog fishery, in particular as water temperatures off the eastern coast of the United States are predicted to continue to increase in response to global climate change. Additionally, this research lends insights into the use of A. islandica growth as a paleoclimate proxy for bottom water temperature. 
    more » « less
  8. Warming in recent decades in the North Atlantic Ocean has been heterogeneous, with locations along the northwestern Atlantic experiencing some of the largest and fastest warming in the last 100 years. This region is important for fisheries but has limited spatial and temporal hydrographic instrumental series extending beyond the past decades, especially along the coastal United States portion of the northwestern Atlantic, thus impacting our understanding of past climatic variability. To provide a longer temporal context for these changes, we constructed a continuous master shell growth chronology spanning the last two centuries and provided geochemical records from the Mid-Atlantic region using the long-lived marine bivalve Arctica islandica. Shells were collected on the outer shelf region off Ocean City, Maryland, in ~ 60 m water depth. This region is sensitive to large-scale North Atlantic Ocean dynamics, including the Atlantic Meridional Overturning Circulation (AMOC) and Gulf Stream eddies. Based on growth histories and shell oxygen isotopes, we provide evidence of hydrographic variability beyond the relatively short instrumental period and evaluate the likely causes for these changes. These data allow us to better characterize recent and past oceanographic changes in the Mid-Atlantic region, synthesize the new results with previously developed paleo-records in the northwestern Atlantic, and provide guidance for the management of fisheries in this region. 
    more » « less